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Abstract 

As a first step to analyze the electromagnetic meson production reactions in the nucleon reso- 
nance region, the parameters of the hadronic interactions of a dynamical coupled-channel model, 
developed in Physics Reports 439, 193 (2007), are determined by fitting the irN scattering data. 
The channels included in the calculations are ttN, rjN and 7r-7riV which has 7rA, pN, and <jN reso- 
nant components. The non-resonant meson-baryon interactions of the model are derived from a set 
of Lagrangians by using a unitary transformation method. One or two bare excited nucleon states 
in each of S, P, D, and F partial waves are included to generate the resonant amplitudes in the 
fits. The parameters of the model are first determined by fitting as much as possible the empirical 
irN elastic scattering amplitudes of SAID up to 2 GeV. We then refine and confirm the resulting 
parameters by directly comparing the predicted differential cross section and target polarization 
asymmetry with the original data of the elastic 7r ± p — > Tr^p and charge-exchange 7r~p — * Tr°n pro- 
cesses. The predicted total cross sections of irN reactions and irN — > nN reactions are also in good 
agreement with the data. Applications of the constructed model in analyzing the electromagnetic 
meson production data as well as the future developments are discussed. 
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* Notice: Authored by Jefferson Science Associates, LLC under U.S. DOE Contract No. DE-AC05- 
06OR23177. The U.S. Government retains a non-exclusive, paid-up, irrevocable, world-wide license to 
publish or reproduce this manuscript for U.S. Government purposes. 

1 



I. INTRODUCTION 



It is now well recognized that a coupled-channel approach is needed to extract the nucleon 
resonance (N*) parameters from the data of ttN and electromagnetic meson production re- 
actions. With the recent experimental developments [H, 0], such a theoretical effort is needed 
to analyze the very extensive data from Jefferson Laboratory (JLab), Mainz, Bonn, GRAAL, 
and Spring-8. To cope with this challenge, a dynamical coupled-channel model (MSL) for 
meson-baryon reactions in the nucleon resonance region has been developed recently [3]. In 
this paper we report a first-stage determination of the parameters of this model by fitting 
the tiN scattering data up to invariant mass W = 2 GeV. 

The details of the MSL model are given in Ref. Here we will only briefly recall its 
essential features. Similar to the earlier works on meson-exchange models j3, 0, @, 0, @, @, 
10|, Q E2, EE H M, M, 0, M, M, 0, H H3, M, H 0, Hf} of pion-nucleon scattering, the 
starting point of the MSL model is a set of Lagrangians describing the interactions between 
mesons (M =7, ir, r\ , p^uj, a, . . .) and baryons (B = N, A, N*, . . .). By applying a unitary 



13, 27 , an effective Hamiltonian is then derived from the considered 



transformation method 
Lagrangian. It can be cast into the following more transparent form 
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Hq + T V + V 22 + h 
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where H = J2a \/ m a + Pa w ith m a denoting the mass of particle a, and 
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(3) 
(4) 



A'* 



MB 



Here c.c. denotes the complex conjugate of the terms on its left-hand-side. In the above 
equations, MB = jN, ttN, rjN, 7rA, pN, aN, represent the considered meson-baryon states. 
The resonance associated with the bare baryon state iV* is induced by the vertex interactions 
T N *^ MB and rjv*-^7r7rAr- Similarly, the bare meson states M* = p, a can develop into reso- 
nances through the vertex interaction h,M*->mr- Note that the masses M%, and m° M , of the 
bare states iV* and M* are the parameters of the model which must be determined by fit- 
ting the tiN and titc scattering data. They differ from the empirically determined resonance 
positions by mass shifts which are due to the coupling of the bare states to the scattering 
states. The term t> 2 2 contains the non-resonant meson-baryon interaction Vmb,m'B' and tttt 
interaction v wn . The non-resonant interactions involving tttcN states are in h^N- All of 
these interactions are energy independent, an important feature of the MSL formulation. 

We note here that the Hamiltonian defined above does not have a ttN <-> N vertex. By 
applying the unitary transformation method, this un-physical process as well as any vertex 
interaction A <-> B + C with a mass relation < itlb + mc are eliminated from the 
considered Hilbert space and their effects are absorbed in the effective interactions v 22 and 
h nn N- This procedure defines the Hamiltonian in terms of physical nucleons and greatly 
simplifies the formulation of a unitary reaction model. In particular, the complications 
due to the nucleon mass and wavefunction renormalizations do not appear in the resulting 
scattering equations. This makes the numerical calculations involving the ttttN channel 
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much more tractable in practice. The details of this approach are discussed in Refs. 13_|, [27 
as well as in the earlier works on wNN interactions [281 ] . 

Starting from the above Hamiltonian, the coupled-channel equations for irN and 'yN 
reactions are then derived by using the standard projection operator technique [291 ]. as given 
explicitly in Ref. (3|. The obtained scattering equations satisfy the two-body (7rN,nN, jN) 
and three-body (irirN) unitarity conditions. The 7rA, pN and <jN resonant components 
of the irnN continuum are generated dynamically by the vertex interaction Ty of Eq. 
Accordingly, the ttttN cuts are treated more rigorously than the commonly used quasi- 
particle formulation within which these resonant channels are treated as simple two-particle 
states with a phenomenological parametrization of their widths. The importance of such a 
dynamical treatment of unstable particle channels was well known in earlier studies of ttN 
scattering 0, 30] and nNN reactions [3ll ]. 

A complete determination of the parameters of the model Hamiltonian defined by Eqs.flTJ- 
(jH) requires good fits to all of the data of nN and 7iV reactions up to invariant mass W < 
about 2 GeV. Obviously, this is a very complex task and can only be accomplished step 
by step. Our strategy is as follows. We need to first determine the parameters associated 
with the hadronic interaction parts of the Hamiltonian. With the fits to nn phase shifts in 
Ref. [32T ]. the nn interactions /i Pj7r7r and and the corresponding bare masses for p and 
a have been determined in an isobar model with v n7T = 0. We next proceed in two stages. 
The first-stage is to determine the ranges of the parameters of the interactions Tn*->mb and 
vmb,m'b'- This will be achieved by fitting the nN scattering data from performing coupled- 
channel calculations which neglect the more complex three-body interaction term h^N- This 
simplification greatly reduces the numerical complexity and the number of parameters to be 
determined in the fits. This first-stage fit will provide the starting parameters to fit both 
the data of nN scattering and nN — ► nnN reactions. In this second-stage, the parameters 
associated with Tn*^,mb and Vmb,m'B' will be refined and the parameters of h nn N are then 
determined. The dynamical coupled-channel calculations for such more extensive fits are 
numerically more complex, as explained in Ref. 0]. 

In this work we report on the results from our first-stage determination of the parameters 
of IV-a/b andv MB ,M'B' of Eqs.©-© with MB, M'B' = nN, rjN, nA, pN, aN. We proceed 
in two steps. We first locate the range of the model parameters by fitting as much as possible 
the empirical nN elastic scattering amplitudes up to W = 2 GeV of SAID (33|. We then 
refine and confirm the resulting parameters by directly comparing our predictions with 
the original nN scattering data. Our procedures are similar to what have been used in 
determining the nucleon-nucleon (NN) potentials [34] from fitting NN scattering data. 

The constructed model can describe well almost all of the empirical nN amplitudes in 
S, P, D, and F partial waves of SAID (33[. We then show that the predicted differential 
cross sections and target polarization asymmetry are in good agreement with the original 
data of elastic n ± p — > n ± p and charge-exchange 7i~p — > 7i°n processes. Furthermore the 
predicted total cross sections of the ttN reactions and ttN — > r]N reactions agree well with 
the data. Thus the constructed model is at least comparable to, if not better than, all of the 



recent nN models [ll|, [12|, [13j, [19|, |20j, [22j, [23|, [2J, [26| . It can be used to perform a first-stage 
extraction of the 7iV — > N* parameters by analyzing the photo- and electro-production of 
single n meson. It has also provided us with a starting point for performing the second-stage 
determination of the model parameters by also fitting the data of 7riV — ► nnN reactions. 
Our efforts in these directions are in progress and will be reported elsewhere. 

In Section II, we recall the coupled-channel equations presented in Ref. j^]. The calcula- 
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FIG. 1: Graphical representation of Eqs.([5])- ([2Tj) . 

tions performed in this work are described in Section III. The fitting procedure is described 
in Section IV and the results are presented in Section V. In Section VI we give a summary 
and discuss future developments. 

II. DYNAMICAL COUPLED-CHANNEL EQUATIONS 



With the simplification that 7T7riV interaction h^rN of Eq. (j3J) is set to zero, the meson- 
baryon (MB) scattering equations derived in Ref. |3| are illustrated in Fig. [TJ Explicitly, 
they are defined by the following equations 



TmB,M'B'(E) = tMB,M'B> M'B> 



(5) 



where MB = irN, r/N, 71 A, pN, aN. The full amplitudes T nN ^ nN (E) can be directly used to 
calculate nN scattering observables. The non-resonant amplitude t MB>M i B i(E) in Eq. (j5J) is 
defined by the coupled-channel equations, 



tMB,M'B'(E) — VmB,M'B'(E) + VmB,M"B"(E) Gm"B"(E) tM"B",M'B'(E) 



M"B" 



with 



VmB,M'B'(E) — VMB,M'B' + Z^MB M'B'{E) ■ 



(6) 



(7) 



Here the interactions vmb,m'b' are derived from the tree-diagrams illustrated in Fig. [2] by 
using a unitary transformation method [l3|, 27]. It is energy independent and free of singu- 
larity. On the other hand, Zj^ B M , B ,(E) is induced by the decays of the unstable particles 
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FIG. 2: Mechanisms for vmb,M'B> of Eq. (TfJ): u s direct s-channel, v u crossed u-channel, v l one- 
particle-exchange t-channel, v c contact interactions. 
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FIG. 3: One-particle-exchange interactions Z^ AwA (E), Z pNnA and Z^ nA of Eq. (JTj). 

(A, p, a) and thus contains moving singularities due to the tttcN cuts, as illustrated in 
Fig 13J Here we note that if the ttttN interaction term h nn N of Eq.fll]) is included, the driving 
term Eq. (j7j) will have an additional term Z^ IB M , B , (E) which involves a 3-3 txtx N amplitude 
t nnNt7T7TN , as given in Ref. [H, and hence is much more difficult to calculate. As explained in 
Section I, we neglect this term in this first-stage fit to the ttN scattering data. 

The second term in the right-hand-side of Eq. (jSJ) is the resonant term defined by 

t MB,M l B'( E )= ^MB^N*{E)[D{E)]ijt N *^ M i B <{E) , (8) 

with 

[D~ l {E)]ij = (E — M° N: )5 hJ - Ey(JE0 , (9) 
where Mjy* is the bare mass of the resonant state N*, and the self-energies are 

^i,j(E) = X! ^ n*^mbG mb{E)Y mb^n*{E) . (10) 

MB * 3 

The dressed vertex interactions in Eq. (jSJ) and Eq. (ITU1) are (defining Tmb-*n* — ^n*^mb) 

FmB->N*{E) = T MB _> N * + tMB,M'B'{E)GM'B'{E)TM'B'->N* , (H) 

M'B' 

Fn*-^Mb(E) = Tn*^mB + X! ^ N*^M'B'GM'B'{E)tM'B' ,Mb{E) . (12) 



A/'B' 



It is useful to mention here that if there is only one N* in the considered partial wave, the 
resonant amplitude (Eq. (JSj)) can be written as 



,R , m _ T MB^N* (E)T N ^M'B'(E) o 

t MB ,M>B>{E)- E _ ER{E)+i r_ J p W 
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with 



E R (E) = M^+Re[Z{E)\, (14) 
F R (E) = -2Im[E( J B)], (15) 



where, 



Fn*->MbGmb(E){ E [SmB,M'B' + tMB,M'B'{E)GM'B'(E)}}r M 'B'^N*(E) . 
MB M'B' 

(16) 

The form Eq. (1131) is similar to the commonly used Breit-Wigner form, but the resonance 
position E R (E) and width T R (E) are determined by the N* — > MB vertex and the non- 
resonant amplitude tMB,M'B'- This is the consequence of the unitarity condition and is 
an important and well known feature of a dynamical approach. Namely, the resonance 
amplitude necessarily includes the non-resonant mechanisms. This feature is consistent 
with the well developed formal reaction theory j29|. Eq. ffTB"]) indicates that it is essential 
to understand the non- resonant mechanisms in extracting the bare vertex functions T^*^mb 
which contain the information for exploring the N* structure. The parameterization used 
for T N * tMB will be explained in Section HTT1 We also note here that the energy dependence of 
E R (E) and T R (E), defined by Eqs f|T11) - ffT5l) . is essential in determining the resonance poles 
in the complex E-plane. 

The meson-baryon propagators Gmb in the above equations are 

E - E M {k) - E B {k) + te 
for the stable particle channels MB = 7rN, r)N, and 

1 

G M B(k,E) = — , r j- r , , (18) 

E - E M {k) - E B (k) - T, MB {k, E) 
for the unstable particle channels MB = it A, pN, aN . The self-energies (36[ in Eq. ( fTBl) are 

m A f 2 M vN (q) \fA,7vN(q)\ 2 



E A (k) J * * [M 2 N (q) + E - E w (k) - [(E N (q) + E^q)) 2 + fc*]i/2 + ze ' 

(19) 



2 , M *M \f P Mi)\' 



Spiv(/C ' E) E p {k) J q dq [MIM + k 2 ] V2 e - E N (k) - [{2EMY + k 2 } 1 / 2 + le ' (20) 

y (k F\ = m<J ( n 2 Jn Mnn ^ \UM<l)\ 2 , 9n 

aN{ ' J 7 9 9 [M^(g) + k 2 } V2 £ - - [(2^(g)) 2 + P]V 2 + *e ' 1 j 

where M nN (q) = E n (q) + E N (q) and M nn (q) = 2E n (q). The vertex function /a^aK?) is 
taken from Ref. [13J], / p 7r7r (g) and L. j7r7r (g) are from the isobar fits [32j to the im phase shifts. 
They are also given explicitly in [3J. 

Here we note that the driving term z[^ B M , B , of Eq. © is also determined by the same 
vertex functions /a.-ttAt (?) ; fp^io) an d fp,nn(l) of Eqs. ffT9l) - fl2Tl) . This consistency is essential 
for the solutions of Eq. (jSJ) to satisfy the unitarity condition. 
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III. CALCULATIONS 



We solve the coupled-channel equations defined by Eqs.(j5|)- (j2ip in the partial- wave repre- 
sentation. The input of these equations are the partial- wave matrix elements of Tn*^mb and 
vmbm'B' of Eqs.©-©, with MB, M'B' = nN,rjN, 7rA, pN, aN, and Z^ BM>B , of Eq. Q 
with MB, M'B' = irA, pN,aN. The calculations of these matrix elements have been given 
explicitly in the appendices of Ref. |3|]. Here we only mention a few points which are needed 
for later discussions. 

In deriving the non-resonant interactions Vmb,m'B' of Eq. ([7]) we consider the tree- 
diagrams (Fig. [2]) generated from a set of Lagrangians with it, rj, a, p, uj, N , and A fields. 
The higher mass mesons, such as a®, a± included in other meson-exchange ttN models, 
such as the Jiilich model are not considered. The employed Lagrangians are ( in the 
convention of Bjorken and Drell 371 ]) 



JttNN 



tNA 



tAA 



j t}NN 



JpNN 



J NNpir 



] ^ N -i>Nl^N-d»^, (22) 



m-n 



h ^^ N -d^, (23) 



rn 

fwAA 
frjNN 



^a^7 5 T a ^a-^4, (24) 
i'Nl^i'Nd' 1 ^ ■ (25) 



9pNn^n[1p ~ ^^°V< 9 1p m • , (26) 



LpNA = -i^^iPaY^T ■ [d ^pl - d v pjj\i}> N + [h.c] , (27) 

Vdp 

LpAA = gpAA^Aah" ~ ^^djfr ■ f^A , (28) 

Lpnir = g p7 rA4>K X 9fj,4> w ] ■ p^ , (29) 
r fnNN 



-gpNN^NlplhTlpN 'P^X^, (30) 



2 

L N Npp = P 9pNN ^NV^TlpN ■ fyX pi. (31) 

8m N 

Lujnn = guNN^Nbn - Tr^—v^d^^N , (32) 

L^ p = - 9 -^e mXv d a p^^ , (33) 

LaNN = gaNN^N^N^a (34) 
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2m n 



(35) 



To solve the coupled-channel equations, Eq. we need to regularize the matrix elements 
of Vmb,m'B'i illustrated in Fig. [2J Here we follow Ref. [13J] in order to use the parameters 
determined in the A (1232) region as the starting parameters in our fits. For the v s and 
v u terms of Fig. [21 we include at each meson-baryon-baryon vertex a form factor of the 
following form 

F(k,A) = [k 2 /[(k 2 + A 2 )] 2 (36) 

with k being the meson momentum. For the meson-meson-meson vertex of v* of Fig. [2j the 
form Eq. (|36|) is also used with k being the momentum of the exchanged meson. For the 
contact term v c , we regularize it by F(k,A)F(k',A'). 

With the non-resonant amplitudes generated from solving Eq. the resonant ampli- 
tude tfy[B m'b' Eq- © then depends on the bare mass M%, and the bare N* — > MB vertex 
functions. As discussed in Ref. jjjj], these bare N* parameters can perhaps be taken from 
a hadron structure calculation which does not include coupling with meson-baryon contin- 
uum states or meson-exchange quark interactions. Unfortunately, such information is not 
available to us. We thus use the following parameterization 



N*,MB(LS) 



(k) 



(27T)3/2 ^ 



c 



N*,MB(LS) 



A 2 



^N*,MB(LS) + ^R) 2 



- (2+L/2) 


" k ' 




_m 7T _ 



(37) 



where L and S are the orbital angular momentum and the total spin of the MB system, 
respectively. The above parameterization accounts for the threshold k L dependence and the 
right power (2 + L/2) such that the integration for calculating the dressed vertex Eq. (iTTj) - 
(|T2|) is finite. Nevertheless as we will discuss in Section [V] this parameterization could be 
too naive. 

The partial- wave quantum numbers for the considered channels are listed in Table [H The 

( F) 

numerical methods for handling the moving singularities due to the tttiN cuts in Z MB M , B , 
(Fig. E]) in solving Eq. (jUJ) are explained in detail in Ref [§]. To get the ttN elastic scattering 
amplitudes, we can use either the method of contour rotation by solving the equations on 
the complex momentum axis k = ke~ ld with 6 > or the Spline-function method developed 
in Refs. [3!§, 39i] and explained in detail in Ref. We perform the calculations using these 
two very different methods and they agree within less than 1%. When Z MB M , B , is neglected, 
Eq. dSJ) can be solved by the standard subtraction method since the resonant propagators, 
Eqs. (Tl8l) . for unstable particle channels irA, pN, and aN are free of singularity on the real 
momentum axis. A code for this simplified case has also been developed to confirm the 
results from using the other two methods. 

The method of contour rotation becomes difficult at high W since the required rotation 
angle 9 is very small. The Spline function method has no such limitation and we can perform 
calculations at W > 1.9 GeV without any difficulty. Typically, 24 and 32 mesh points 
are needed to get convergent solutions of the coupled-channel integral equation ([6]). Such 
mesh points are also needed to get stable integrations in evaluating the dressed resonance 
quantities Eqs. (fT0D-(fT2l). 
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(LS) of the considered partial waves 





ttN 


rjN 


7rA 


criV 


pN 




(0,|) 


(0,|) 


(2,1) 


(1,1) 


(0,1), (2, 


2> 


5*31 


(o,|) 


- 


(2,1) 


- 


(0,|), (2, 


-) 

2> 


Pu 


(1,1) 


(1,1) 


(1,1) 


(o,|) 


(1,|), (1, 


^) 
2> 


Pl3 


(1,1) 


(1,1) 


(i,|),(3,D 


(2,i) 


(i, i),(i,D, 


(3,|) 


P S 1 


(1,1) 


— 


(1,1) 


— 


(i,|), a, 


2> 


-P33 


(1,1) 


- 


(1,|),(3,|) 


- 


(1, i),(i, I), 


(3,|) 


D 13 


(2,|) 


(2,*) 


(0,|),(2,|) 


(1,1) 


(2,|), (0,|), 


(4,1) 


D15 


(2,|) 


(2,|) 


(2,|), (4,|) 


(M) 


(2, 2), (2, 2)' 


(4,1) 


D 33 


(2,|) 




(0,|),(2,|) 




(2,|), (0,|), 


(2,1) 


D 35 


(2,|) 




(2,|), (4,|) 




(2, 2), (2, 2)' 


(4,1) 


Fl5 


(3,|) 


(3,|) 


(1,|),(3,|) 


(2,|) 


(3, ^), (1, |), 


(3,1) 


F17 


(M) 


(3,|) 


(3, 2), (5, 2) 


(4,|) 


(3, 2), (3, 2)' 


(5,|) 


F35 


(3,1) 




(1,|),(3,|) 




(3, 2), (1, 2)' 


(3,1) 


F37 


(3,|) 




(3,|),(5,|) 




(3, 2), (3, 2), 


(5,1) 



TABLE I: The orbital angular momentum (L) and total spin (S')of the partial waves included in 
solving the coupled channel Equation ([6|). 



IV. FITTING PROCEDURE 

With the specifications given in Section III, the parameters associated with Z\ rg AI , B , of 
Eq. are completely determined from fitting the tttt phase shifts in Refs. [13| and [32| . Thus 
the considered model has the following parameters: (a) the coupling constants associated 
with the Lagrangians listed in Eqs. ( |22l) -( l35l) . (b) the cutoff A for each vertex of vmb,m'B' 
(Fig. [2]), (c) the coupling strength Cn\mb{ls) and range kn and ^n*,mb{ls) of the bare 
N* — > MB vertex Eq. ( |37|) . and (d) the bare mass Mjy, of each A^* state. We determine 
these by fitting the ttN scattering data. 

Our fitting procedure is as follows. We first perform fits to the ttN scattering data up to 
about 1.4 GeV and including only one bare state, the A (1232) resonance. In these fits, the 
starting coupling constant parameters of vmb,m'B' are taken from the previous studies of 
ttN and A^ A" scattering, which are also given in Ref . 0] . Except the 71 NN coupling constant 
fnNN all coupling constants and the cutoff parameters are allowed to vary in the x 2_n t to 
the 71N data. The coupled-channel effects can shift the coupling constants greatly from their 
starting values. We try to minimize these shifts by allowing the cutoff parameters to vary 
in a very wide range 500 MeV < A < 2000 MeV. Some signs of coupling constants, which 



could not be fixed by the previous works 40], are also allowed to change. We then use the 



parameters from these fits at low energies as the starting ones to fit the amplitudes up to 
2 GeV by also adjusting the resonance parameters, M%,, Cn*,mb(ls), kn and A.jsr*,MB(LS)- 
Here we need to specify the number of bare N* states in each partial wave. The simplest 
approach is to assume that each of 3-star and 4-star resonances listed by the Particle Data 
Group [35[ is generated from a bare N* state of the model Hamiltonian Eq. ([1]). However, 
this choice is perhaps not well justified since the situation of the higher mass N*'s is not so 
clear. 



9 



We thus start the fits including only the bare states which generate the lowest and well- 
established N* resonance in each partial wave. The second higher mass bare state is then 
included when a good fit can not be achieved. We also impose the condition that if the 
resulting M%, is too high > 2.5 GeV, we remove such a bare state in the fit. This is due 
to the consideration that the interactions due to such a heavy bare N* state could be just 
the separable representation of some non-resonant mechanisms which should be included 
in vmb,m'b'- In some partial waves the quality of the fits is not very sensitive to the N* 
couplings to 7rA, pN, and aN. But the freedom of varying these coupling parameters is 
needed to achieve good fits. 

It is rather difficult to fit all partial waves simultaneously because the number of resonance 
parameters to be determined is very large. We proceed as follows. We first fit only 3 or 
4 partial waves which have well established resonant states, and whose amplitudes have an 
involved energy dependence. These are the Sn, Pn, S31 and P 33 partial waves. These fits are 
aimed at identifying the possible ranges of the parameters associated with Vmb,m'b'- This 
step is most difficult and time consuming. We then gradually extend the fits to include more 
partial waves. For some cases, the fits can be reached easily by simply adjusting the bare 
iV* parameters. But it often requires some adjustments of the non-resonance parameters to 
obtain new fits. This procedure has to be repeated many times to explore the parameter 
space as much as we can. We carry out this very involved numerical task by using the fitting 
code MINUIT and the parallel computation facilities at NERSC in US and the Barcelona 
Super computing Center in Spain. 

The most uncertain part of the fitting is to handle the large number of parameters asso- 
ciated with the bare N* states. Here the use of the empirical partial-wave amplitudes from 
SAID is an essential step in the fit. It allows us to locate the ranges of the N* parameters 
partial- wave by partial- wave for a given set of the parameters for the non- resonant vmb,m'b,- 
Even with this, the information is far from complete for pinning down the N* parameters. 
Perhaps the N* parameters associated with the ttN state are reasonably well determined in 
this fit to the irN scattering data. The parameters associated with rjN, it A, pN and aN 
can only be better determined by also fitting to the data of ttN — > rjN and ttN —>■ irnN 
reactions. This will be pursued in our second-stage calculations, as discussed in section I. 

It is useful to note here that the leading-order effect due to Z^ E > of the meson-baryon 
interaction Eq. (j7|) on nN elastic scattering is 

Svkn^n — v 7T n,mbGmb{E)Z^ b ^ M i B iGm'b , {E)v M i B i 7tN . (38) 

MB,M'B'=irA,pN,<rN 

We have found by explicit numerical calculations that Sv^n,ttN is much weaker than v^n^n 

( F) 

and hence the coupled channel effects due to Z MB M , B , on ttN elastic scattering amplitude 
are weak. One example obtained from our model is shown in Table [Til Thus we first 
perform the fits without including term to speed up the computation. We then refine 
the parameters by including this term in the fits. 

V. RESULTS 

As mentioned in section I, we first locate the range of the parameters by fitting the 
empirical ttN scattering amplitude of SAID [33]. We then check and refine the resulting 
parameters by directly comparing our predictions with the original irN scattering data. 
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Im[t nN>nN (ZW = 0)] 


Sn 


-0.00481 


-0.00557 




0.0841 


0.0827 


Pn 


0.0937 


0.103 




0.636 


0.640 


Pl3 


0.169 


0.181 




0.275 


0.275 


D l3 


0.202 


0.194 




0.299 


0.309 


D 15 


0.117 


0.116 




0.0179 


0.0179 


Fi 5 


0.290 


0.291 




0.157 


0.155 


F l7 


0.0360 


0.0359 




0.00293 


0.00289 


S31 


-0.433 


-0.437 




0.496 


0.504 


P31 


-0.253 


-0.230 




0.434 


0.448 


1 ) 

-T33 


U.U0U0 


U.U3U6 




U.51U 


U.45 1 


£>33 


-0.00504 


-0.0135 




0.106 


0.104 


D 35 


0.0551 


0.0551 




0.0540 


0.0537 


F35 


-0.0214 


-0.0229 




0.0259 


0.0283 


F37 


0.0625 


0.0626 




0.00502 


0.00512 



TABLE II: The effect of Z MB M , B , on the ttN scattering amplitudes t w N t7r N from solving Eq. ([6]) 
at W = 1.7 GeV. The normalization is t^x^N = (e 2tS7rN — l)/(2i), where S^n is the ttN scattering 
phase shift which could be complex at energies above the ir production threshold. 



Our fits to the empirical amplitudes of SAID [33| are given in Figs. HH5] and Figs. [M3 for 
the T = 1/2 and T = 3/2 partial waves, respectively. The resulting parameters are presented 
in Appendix I. The parameters associated with the non- resonant interactions, vmb,m'b' with 
MB,M'B' = 7TN,r]N, irA, pN, aN, are given in Table HTT1 for the coupling constants of the 
starting Lagrangian Eqs. (l22l) - fl3"o1) and Table HVl for the cutoffs of the form factors defined 
by Eq. fl3"oT) . The resulting bare A^* parameters are listed in Tables IVTlVIII 

From Figs. HE, one can see that the empirical irN amplitudes can be fitted very well. 
The most significant discrepancies are in the imaginary part of S31 in FigJTl The agreement 
is also poor for the Fn in FigHJl5] and -D35 in FigsJMZl but there are rather large errors in 
the data. Our parameters are therefore checked by directly comparing our predictions with 
the data of differential cross sections da/dfl and target polarization asymmetry P of elastic 
n ± p — > ir^p and charge-exchange n~p — > 7r°n processes. Our results (solid red curves) are 
shown in FigsJ8]fT2"l Clearly, our model is rather consistent with the available data, and are 
close to the results (dashed blue curves) calculated from the SAID's amplitudes. Thus our 
model is justified despite the differences with the SAID's amplitudes seen in FigSJH 

It will be important to further refine our parameters by fitting the data of other ttN 
scattering observables, such as the recoil polarization and double polarization. Hopefully, 
such data can be obtained from the new hadron facilities at JPARC in Japan. 

Our model is further checked by examining our predictions of the total cross sections a tot 
which can be calculated from the forward elastic scattering amplitudes by using the optical 
theorem. The total elastic scattering cross sections o el can be calculated from the predicted 
partial wave amplitudes. With the normalization < k\k' >= 5(k — k') used in Ref. [s], we 
have 

a e \W) = ct t W ( 39 ) 

T=l/2,3/2 
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FIG. 4: Real parts of the calculated ttN partial wave amplitudes (Eq. ([5])) of isospin T 
compared with the energy independent solutions of Ref . [33 ] . 
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FIG. 5: Imaginary parts of the calculated ttN partial wave amplitudes (Eq. ([5])) of isospin T 
are compared with the energy independent solutions of Ref. 33] . 
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FIG. 6: Real parts of the calculated ttN partial wave amplitudes (Eq. ([5])) of isospin T = 3/2 are 
compared with the energy independent solutions of Ref . [33 ] . 



with 



K JLS Z 



(k,k,W)f 



(40) 



where p n N(W) = 7ikE n (k)E N (k)/W with k determined by W = E^k) + E N (k) and the 
amplitude T^,,^ LS , N dk,k',W) is the partial- wave solution of Eq. Similarly, the 
total 7rN — > rjN cross sections can be calculated from 



tot 
T TvN^r]N 



JLS Z 



A- 2 



(41) 



where p vN (W) = irk'E v (k')E N (k')/W with k' determined by W = E v (k') + E N (k'). We 
can also calculate the contribution from each of the unstable channels, ttA, pN, and crN, 
to the total tcN — > ttttN cross sections. For example, we have for the ttN — > 7rA — > irnN 
contribution in the center of mass frame 



W — VTLt, 



dM. 



M, 



N 



nN 



r„ A {k,E)/{2ir) 



E A (k) \W - E«(h) - E A (k) - S^ A (fc, E) 



;(T«N^*A(k,W) (42) 



where k is defined by M nN = E n (k) + E N (k), E nN (k) = [M% N + /c 2 ] 1/2 , S 7rA (A:, E) is defined 
in Eq.(HSD, F^ A (k,E) = -2Im(E nA (k, E)), and 



2 J + 1 

&irN-*vA{k, W) = 47rp 7rA r(A; )p 7rA (A;) (25^ + 1)(25 + 1) ^ A ( L ' S ') 



L'S'LS.J 



n(ls) (k,k ;W)\ 2 
(43) 



where k is defined by W = E n (k ) + E N (k Q ) and p a b(k) = 7tkE a (k)E b (k)/W . The am- 
plitude T[, s ,,£s ls(ttN) ^o! WO is the partial-wave solution of Eq.©. The corresponding 
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FIG. 7: Imaginary parts of the calculated irN partial wave amplitudes (Eq. ©) of isospin T 
are compared with the energy independent solutions of Ref. [33]. 
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FIG. 8: Differential cross section for several different center of mass energies. Solid red curve 
corresponds to our model while blue dashed lines correspond to the SP06 solution of SAID 331 ] . 
All data have been obtained through the SAID online applications. Ref. 331 ] . 
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FIG. 9: Differential cross section for several different center of mass energies. Similar description 
as Fig. [HJ All data have been obtained through the SAID online applications. Ref. [331 ]. 
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FIG. 10: Target polarization asymmetry, P, for several different center of mass energies. Similar 
description as Fig. [51 All data have been obtained through the SAID online applications. Ref. [331 ] . 
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FIG. 11: Target polarization asymmetry, P, for several different center of mass energies. Descrip 
tion as in Fig. [5J All data have been obtained through the SAID online applications. Ref. 331 ] 



expressions for the unstable channels pN and crN can be obtained from Eqs. (I4"2l - (l4"3"l) by 
changing the channel labels. 

The predicted a tot (solid curves) along with the resulting total elastic scattering cross 
sections a el compared with the data of 7T + p reaction are shown in Fig. [13j Clearly, the 
model can account for the data very well within the experimental errors. Here only the 
T = 3/2 partial waves are relevant. Equally good agreement with the data for ir~p reaction 
are shown in the left side of Fig. [TH In the right side, we show how the contributions from 
each channel add up to get the total cross sections. The comparison of the contribution 
from nN channel with the data is shown in Fig. [Jj3 It is possible to improve the fit to this 
data by adjusting N* — > rjN parameters. But this can be done correctly only when the 
differential cross section data of tcN —>■ r]N are included in the fit. This is beyond the scope 
of this work and will be pursed in our second-stage calculations. 

The contributions from 7rA, pN and aN intermediate states to the ir~p — > ttttN total 
cross sections calculated from our model can be seen in the right side of Fig. [TH These 
predictions remain to be verified by the future experiments. The existing ttN — > nirN data 
are not sufficient for extracting model independently the contributions from each unstable 
channel. 

The results shown in Figs. [TBlfTol indicate that our parameters are consistent with the 
total cross section data. 

We now discuss the parameters presented in Appendix A. It is rather difficult to compare 
the resulting non- resonant coupling constants listed in Table IIHI with the values from other 
works, since the coupling strengths are also determined by the cutoff parameters listed in 
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FIG. 12: Target polarization asymmetry, P, for several different center of mass energies. Descrip 
tion as in Fig. [HJ All data have been obtained through the SAID online applications. Ref. 33|] 



Table IIVI Perhaps it is possible to narrow their differences by using a different parameter- 
ization of the form factors. However, the fit is a rather time consuming process and hence 
no attempt is made in this work to try other forms of form factors. 

In Table |Vl we see that all of the bare masses are higher than the PDG's resonance 
positions. This can be understood from the expression Eq. (fH|) for the partial waves with 
only one N* since one finds in general that Re[Y>(E)] < 0. For the Sn, Pn, P33 and D ri 
partial waves, two bare N* states are mixed by their interactions, as can be seen in Eq. (jTUl) . 
Thus the relation between their bare masses and the resonance positions identified by PDG 
is much more complex. 

As we mentioned above, the fit to ttN elastic scattering is can not determine well the 
bare A^* — > 7rA, pN,aN parameters. Thus the results for these unstable particle channels 
listed in Tables IIimVHI must be refined by fitting the tcN — ► ttttN data. 



VI. SUMMARY AND FUTURE DEVELOPMENTS 

Within the formulation developed in Ref. 0] , we have constructed a dynamical coupled- 
channel model of ttN scattering by fitting the nN scattering data. The parameters of the 
model are first determined by fitting as much as possible the empirical irN elastic scattering 
amplitudes of SAID up to 2 GeV. We then refine and confirm the resulting parameters by 
directly comparing the predicted differential cross section and target polarization asymmetry 



with the original data of the elastic n p 



7T 



p and charge-exchange it p — > n n processes. 
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FIG. 13: The predicted total cross sections of the ir + p — > X (solid curve)and n + p — * ir + p (dashed 
curve) reactions are compared with the data. Squares and triangles are the corresponding data 



from Ref. 35] 



The predicted total cross sections of irN reactions and are also in good agreement with the 
data. The model thus can be used as a starting point for analyzing the very extensive data 
of electromagnetic it production reactions. 

The predicted total cross sections of ttN — > r]N reactions are also in fair agreement with 
the data. However, the parameters associated with the r)N channel need to be refined to 
also fit the differential cross section data of ttN — > t]N before the model can be used to 
analyze the data of electromagnetic rj production reactions. 

The main shortcoming of this work is that the ttttN interaction term h^N of Eq. (j4j) is 
not included in the calculations. As derived in Ref. 0], the effects due to this interaction 
can be included by adding a term Z$b m'B'(E), which contains the ttttN — ► tttyN scattering 
amplitude, to the driving term Vmb,m'b'{E) of Eq.flB]). Our effort in this direction is in 
progress along with the development of a more complete determination of the parameters 
of the model by fitting both the data of irN elastic scattering and itN — > ttttN reactions. 
This is also essential to pin down the parameters of the interactions associated with the 
7i A, pN and aN states. Only when this second-stage is completed, we then can perform 
dynamical coupled- channel analysis of the very extensive and complex data of photo- and 
electro-production of two pions. This is an essential step to probe the W > about 1.7 GeV 
resonance region where the information on N* is very limited and uncertain. 

Finally, a necessary next step is to extract the resonance poles and the associated residues 
from the predicted ttN amplitudes. This is being pursued and will be published else- 
where (4r]| . 
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APPENDIX A: PARAMETERS FROM THE FITS 



Parameter 




SL Model 






0.08 


0.08 




m a (MeV) 


500.1 






fnNA 


2.2061 


2.0490 




frjNN 


3.8892 






QpNN 


8.7214 


6.1994 




Kp 


2.654 


1.8250 




9wNN 


8.0997 


10.5 






1.0200 


0.0 




9crNN 


6.8147 






9p7T7V 


4. 


6.1994 




/vrAA 


1.0000 






fpNA 


7.516 








2.353 






dump 


6.955 






9pAA 


3.3016 






kpAA 


2.0000 







TABLE III: The parameters associated with the Lagrangians Eqs.(22)-(35). The results are from 
fitting the empirical ttN partial- wave amplitudes 33J of a given total isospin T = 1/2 or 3/2. The 
parameters from the SL model of Ref. 13J are also listed. 
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Parameter 


(MeV) 


SL model (MeV) 




A-nNN 


809.05 


642.18 




AttATA 


829.17 


648.18 




A-pNN 


1086.7 


1229.1 




Ap7T7r 


1093.2 


1229.1 




A-wNN 


1523.18 


— 




A-rjNN 


623.56 


— 




A-aNN 


781.16 


— 




A 

A pNA 


1200.0 






A,rAA 


600.00 








1200.0 








600.00 






ApAA 


600.00 







TABLE IV: Cut-offs of the form factors, Eq. (|36p . of the non-resonant interaction vmb,M'B'- The 
results are from fitting the empirical irN partial-wave am plit udes 331 ] of a given total isospin 
T = 1/2 or 3/2. The parameters from the SL model of Ref. 13| are also listed. 



Ltj 


PDG 


s Mass( MeV) 


Mi (MeV) 


M 2 (MeV) 


Sn 


1535; 


1655 


1800. 


1880. 


S 3 i 


1630 




1850. 




Pu 


1440; 


1710 


1763 


2037 


Pl3 


1720 




1711 




Psi 


1910 




1900.3 




P33 


1232; 


1600 


1391 


1602. 


D 13 


1520; 


1700 


1899.1 


1988. 


D 15 


1675 




1898 




D33 


1700 




1976 




D 35 


1960 








F15 


1685 




2187 




F35 


1890 




2162 




F37 


1930 




2137.8 





TABLE V: The masses of the nucleon excited states included in the fits, (second and third 
columns). The first column contains the masses of the nucleon resonances given by PDG [35]. 
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7riV 


nN 


7rA 


aN 


pN 


Sn (1) 
uil K^J 


7.0488 


9.1000 


-1.8526 




—2.7945 


2.0280 


.02736 




Sn (2) 


9.8244 


60000 


.04470 




1.1394 


-9.5179 


-3.0144 




Sn 


5 275002 




-6.17463 






-4.2989 


5.63817 






3.91172 


2 62103 


-9.90545 




—7.1617 


-5.1570 


3.45590 




Pn (2) 
± 11 l^v 


9.9978 


3.6611 


-6.9517 




8 62949 


-2.9550 


-0.9448 




Pii 
1 13 


3.2702 


-.99924 


-9.9888 


-5.0384 


1.0147 


-.00343 


1.9999 


-.08142 


- 1 31 


6 80277 




2.11764 






9.91459 


0.15340 






1 31883 




2.03713 


9.53769 




-.3175 


1.0358 


0.76619 


P« (2) 
1 33 


1.3125 




1.0783 


1.52438 




2.0118 


-1.2490 


0.37930 


-^13 V x 7 


.44527 


-.0174 


-1.9505 


.97755 


- 481855 


1.1325 


-.31396 


.17900 


Die f2l 

-^13 


46477 


35700 


9.9191 


3.8752 


—5 4QQ4 


.28916 


9.6284 


-.14089 


015 


.31191 


-.09594 


4.7920 


.01988 


-.45517 


-.17888 


1.248 


-.10105 


033 


.9446 




3.9993 


3.9965 




.16237 


3.948 


-.85580 


05 


.06223 


0.0000 


1.0395 


.00454 


1.5269 


-1.0353 


1.6065 


-.0258 


035 


.173934 




-2.96090 


-1.09339 




-.07581 


8.0339 


-.06114 


037 


0.25378 




-0.3156 


-0.0226 




.100 


.100 


.100 



TABLE VI: The coupling constants C N *, JT LS;MB of Eq. ([37]) with MB = irN, rjN, vrA, aN, pN 
for each of the resonances. When there are more than one value for irA and pN channels, they 
correspond to the possible quantum numbers (LS) listed in Table 2. 





irN 


rjN 


7rA 


aN 


pN 


Sn (i) 


1676.4 


598.97 


554.04 




801.03 


1999.8 


1893.6 




511 (2) 


533.48 


500.02 


1999.1 




1849.5 


796.83 


500.00 




5*31 


2000.00 




500.00 






500.031 


500.00 




Pn (1) 


1203.62 


1654.85 


729.0 




1793.0 


621.998 


1698.90 




Pn (2) 


646.86 


897.84 


501.26 




1161.20 


500.06 


922.280 




Pl3 


1374.0 


500.23 


500.00 


500.770 


640.50 


500.00 


500.10 


1645.2 


P31 


828.765 




1999.9 






1998.8 


2000.6 




P33 (1) 


880.715 




507.29 


501.73 




606.78 


1043.4 


528.37 


P33 (2) 


746.205 




846.37 


780.96 




584.98 


500.240 


1369.7 


013 (1) 


1658. 


1918.2 


976.36 


1034.5 


1315.8 


599.79 


1615.1 


1499.50 


013 (2) 


1094.0 


678.41 


1960.0 


660.02 


1317.0 


550.14 


597.57 


1408.7 


015 


1584.7 


1554.0 


500.77 


820.17 


507.07 


735.40 


749.41 


937.53 


033 


806.005 




1359.38 


608.090 




1514.98 


1998.99 


956.61 


015 


1641.6 


655.87 


1899.5 


522.68 


500.93 


500.76 


500.0 


1060.9 


035 


1035.28 




1227.999 


586.79 




1514.3 


593.84 


1506.0 


037 


1049.04 




1180.2 


1031.81 




600.02 


600.00 


600.02 



TABLE VII: The range parameter An*,jtls-,mb (in unit of (MeV/c)) of Eq. ([37]) with MB = 
irN, rjN, ttA, aN, pN for each of the resonances. When there are more than one value for rrA and 
pN channels, they correspond to the possible quantum numbers (LS) listed in Table 2. 
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